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Abstract 

This paper presents an investigation into spreading dynamics and dynamic contact angle of Ti0 2 -deionized water 
nanofluids. Two mechanisms of energy dissipation, (1) contact line friction and (2) wedge film viscosity, govern the 
dynamics of contact line motion. The primary stage of spreading has the contact line friction as the dominant 
dissipative mechanism. At the secondary stage of spreading, the wedge film viscosity is the dominant dissipative 
mechanism. A theoretical model based on combination of molecular kinetic theory and hydrodynamic theory 
which incorporates non-Newtonian viscosity of solutions is used. The model agreement with experimental data is 
reasonable. Complex interparticle interactions, local pinning of the contact line, and variations in solid-liquid 
interfacial tension are attributed to errors. 

Keywords: Dynamic contact angle, Hydrodynamic theory, Molecular kinetic theory, Nanofluids, Nanoparticles, 
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Background 

Industrial operations such as spin coating, painting, and 
lubrication are based on spreading of fluids over solid 
surfaces. The fluid may be simple [1-3] or particulate 
such as paint, ink, or dye [4]. For many years, capillary 
flow of simple fluids has received considerable attention, 
and physics of capillary action is known for a long time 
[5-9]. In addition, capillary flow of micellar surfactant 
solutions which contain monodisperse and naturally sta- 
bilized nanoparticles has been studied [10-14], However, 
the same study on liquids laden with metallic and oxide 
nanoparticles such as silver, copper, zinc oxide, and ti- 
tanium oxide is scarce. These fluid suspensions are 
termed as nanofluids after the seminal work by Choi and 
Eastman [15]. The application of nanofluids is coined 
with enhanced heat transfer performance compared with 
their base fluids. They are proposed for applications in 
cooling of electronic devices, ventilation and air condi- 
tioning, and biomedical applications [14,16-24]. 

It is known that out-of-equilibrium interfacial energy 
(<r(cos 6° - cos 0)) provides free energy of capillary flow 
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where a is the liquid-air surface tension and 6° and 6 are 
the equilibrium and dynamic contact angles, respect- 
ively. During capillary flow, the free energy is dissipated 
by two mechanisms [5]: (1) contact line friction (TZ/) 
which occurs in proximity of three-phase contact line 
(solid-liquid-air). The friction at the three-phase con- 
tact line is due to intermolecular interactions between 
solid molecules and liquid molecules. (2) Wedge film 
viscosity {TZ W ) which occurs in the wedge film region 
behind the three-phase contact line. Lubricating and 
rolling flow patterns in the wedge film region result in 
the dissipation of the free energy. For each mechanism 
of energy dissipation, a theory is developed: (1) molecu- 
lar kinetic theory (MKT) [25,26] models the contact line 
friction, and (2) hydrodynamic theory (HDT) [27,28] 
models the wedge film viscosity. For partial wetting sys- 
tems (6° > 10°), it is assumed that both dissipative mech- 
anisms coexist and models that combine MKT and HDT 
are developed by Petrov [29] and De Ruijter [30]. In 
Petrov's model, it is assumed that the equilibrium con- 
tact angle 6° is not constant and its change is described 
by MKT. In De Ruijter s model, it is assumed that 6° is 
constant and the dissipation functions are added to form 
the total dissipation function, TZ tot = + TZ W . These 
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models are developed for Newtonian fluids and show 
generally good agreement with experimental data [31]. 

This paper presents an investigation into spreading 
dynamics and dynamic contact angle of Ti0 2 -deionized 
(DI) water nanofluids. Metal oxide Ti0 2 nanoparticle 
was chosen for its ease of access and popularity in en- 
hanced heat removal applications. Various nanoparticle 
volume concentrations ranging from 0.05% to 2% were 
used. The denser solutions exhibit non-Newtonian vis- 
cosity at shear rate ranges that are common to capillary 
flow. To model experimental data a theoretical model 
based on combination of MKT and HDT similar to De 
Ruijter's model is used. The non-Newtonian viscosity of 
the solutions is incorporated in the model. 

Methods 

Preparation of nanofluids 

The solutions were prepared by dispersing 15 nm 
Ti0 2 nanoparticles (anatase, 99%, Nanostructured and 
Amorphous Materials Inc., Houston, TX, USA) in DI water. 
Oleic acid is reported to stabilize Ti0 2 nanoparticles in DI 
water [20] and was added to the mixture at 0.01vol.% con- 
centration. The solution was stirred for 8 h followed by 100 
min sonication (Sonicator 3000, 20 kHz and 80 kW, 
MISONIX, Farmingdale, NY, USA). Temperature of the 
solution was maintained at 25°C during the sonication 
process. Clustering and morphology of nanoparticles are 
important factors in nanofluid spreading capability. We 
used transmission electron microscopy (TEM) nanographs 
of Ti0 2 nanoparticles to examine these factors. An aliquot 
of dilute solution was dropped and dried on a carbon- 
coated copper grid. TEM images were then taken immedi- 
ately. Figure 1 shows that the solution contains irregular 
particle clusters in addition to monodispersed particles. 
The sizes of the single particles were found to be close to 
15 nm as specified by the supplier. The morphology of the 




Figure 1 TEM nanographs of 15 nm Ti0 2 nanoparticles. 
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monodispersed particles is spherical. Sonication of the 
nanofluid solution and addition of surfactant molecules is 
critical to break down the particle agglomerations and 
stabilize particle dispersion. The effective nanoparticle 
size was 260 nm measured with a particle size analyzer 
(Brookhaven Instruments Corporation, Holtsville, NY, 
USA). Adsorption of oleic acid surfactant molecules to the 
surface of Ti0 2 particles and dissociation of proton from 
carboxylic acid head groups result in net negative charges 
on the surface of particles and thus formation of electric 
double layer around them. Thick electric double layers 
cause the deviation of particle-particle interactions from 
hard-sphere interactions. The (Debye) length in nanometer 
of an electric double layer of 1:1 electrolyte in water at 25°C 
can be approximated by 0.3/ \[M (where M is the molar 
concentration). For 0.01 vol.% concentration of oleic acid in 
water (which is 3.15 x 10" 4 molar), the Debye length is 
estimated to be about 16.9 nm. Such a small increase in the 
effective diameter of particles allows for an assumption of 
hard-sphere interactions between particles in the solution 
which is an important assumption in using Krieger's 
formula [32]. All other experimental measurements were 
carried out at 25°C. 

Measurement of viscosity 

Viscosity of the solutions was measured using a control- 
lable low shear rate concentric cylinders rheometer 
(Contraves, Low Shear 40, Zurich, Switzerland). The 
viscosity was measured at shear rates ranging from 0 to 
50 s 1 . This range corresponds to the shear rates that 
are common to capillary flow. 

Measurement of surface tension 

Surface tension of the solutions was measured by pen- 
dant droplet method using FTA200 system (First Ten 
Angstroms, Inc., Portsmouth, VA, USA). To form the 
pendant droplets, the solutions were pumped out of a 
syringe system at a very low rate, namely 1 ul/s, to 
minimize inertia effects. To minimize errors due to 
evaporation, surface tension was measured right after 
the pendant droplet reached its maximum volume, 
namely 10 ul for the dense solutions. 

Measurement of dynamic contact angle 

Dynamic contact angle of the solutions was measured 
using the FTA200 system. A droplet of solution was gen- 
erated at a very low rate (1 ul/s) and detached from the 
syringe needle tip as soon as it touched the borosilicate 
glass slide. A video was captured while the droplet was 
spreading over the glass slide from initial contact to 
equilibrium position (see Figure 2: the time frame 
elapses between (a) to (b) and (b) to (c) are 5 and 100 s, 
respectively). The consecutive photographs were used to 
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measure the contact angles. The spatial resolution was 
estimated to be about 50 \im on the basis of the focused 
area and camera pixel size. The standard deviation for 
contact angle measurements was less than 1°. The tem- 
poral resolution was estimated based on the frame speed 
of the CCD camera as 30 fps. For each concentration, 
three experiments were performed and average was 
taken. 

Theory 

Empirical analysis of viscosity 

From Figure 3, it is obvious that 0.5%, 1%, and 2% solu- 
tions exhibit shear thinning viscosity at shear rates 
below 20 s" 1 . At higher shear rates, Newtonian behavior 
was observed for all solutions. For dilute solutions, 0.1 
vol.% and 0.05 vol.%, a weak shear thinning behavior 
was also observed at very low shear rates [19]. 

A power-law equation is used to model the shear rate 
and nanoparticle concentration dependent viscosity: 



F(cf>)Ky n 



(1) 



where rj b is the viscosity of DI water equal to 0.927 
mPa s, F((p) is a function of nanoparticle volume 



concentration (0), Ky n ~ l is an indicator of shear thin- 
ning viscosity with K as the proportionality factor, and 
n as the power-law index. F((p) is calculated using 
Krieger's formula [32]: 



1- 



4> 



-2.5<P„ 



(2) 



where 0 max is the fluidity limit that is empirically equal 
to 0.68 for hard spherical particles. In Equation 1, n 
and K are empirical constants which are obtained by 
fitting this equation to the experimental data shown in 
Figure 3. Table 1 shows the values of K and n for vari- 
ous nanoparticle volume concentrations. It is obvious 
that higher nanoparticle concentration results in a lar- 
ger non-Newtonian behavior. Figure 3 also shows that 
the power-law Equation 1 is in good agreement with 
the experimental data. 

Molecular kinetic theory 

Schematic of a spreading droplet of radius r and contact 
angle 0 that is inspired by De Gennes [5] and Blake [26] 
is depicted in Figure 4. Based on MKT [26], the rate of 
displacement of the three-phase contact line over ad- 
sorption sites on solid surface, U, is equal to the net 
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Table 1 Power-law viscosity, surface tension, and equilibrium contact angle of Ti0 2 -Dl water solutions 



Ti0 2 volume concentration (cp) Power-law index (n) Proportionality factor {K) Surface tension (a [N/m]) Equilibrium contact angle (0°) 



2% 


0.04 


2,932 


0.0543 


51.7 


1% 


0.18 


432 


0.0606 


47.5 


0.5% 


0.76 


5 


0.0612 


46.7 


0.1% 


0.89 


2 


0.0623 


45.7 


0.05% 


0.92 


1 


0.0632 


44.5 



frequency of molecular movements, K w (K w = IC - IC, 
where IC is the frequency of forward motion and IC is 
the frequency of backward motion), multiplied by aver- 
age distance between the adsorption sites, A: 



U = AKy 



(3) 



The equilibrium frequency of the three-phase contact 
line motion {K° w ) is obtained from Eyring's theory of ab- 
solute reaction rates [26]: 



IC 



exp 



-AG 



w 



Nk B T, 
At equilibrium K° w 



(4) 



IC 



where k B ,h^N, and T are the Boltzmann constant, 
Planck constant, Avogadro's number, and absolute 
temperature, respectively. In this equation, AG* W is the 
equilibrium free energy of capillary flow. An imbalance 
of the three interfacial tensions near the three-phase 
contact line, solid-liquid (<r s i), solid- vapor ((X sa ), and li- 
quid-vapor (<x), results in the out-of-equilibrium inter- 
facial energy (cr(cos 6° - cos 6)) which changes the total 
free energy of capillary flow. The frequency of the three- 
phase contact line motion in forward direction (+) and 
backward direction (-) is [26]: 



hJ eXP \Nk B T 



AG* W ^ <r( cos#°- cos#) 



2nk B T 



(5) 



where n is the number of adsorption sites per unit area 
on solid surface. The net frequency of contact line mo- 
tion is then as follows [26]: 



K w = K + -IC 



^ exp ( ) 



2K 

T)^Nk B T 



^(cos^-costf)" 
k 2nk B T 



(6) 



For small arguments of sinh, Equations 3 and 6 result 
in linear MKT [31]: 



a(cos0°-cos0) = (U 



(7) 



where C = nk B T / 'K° W X is in units of Pa s and is termed 
as the coefficient of friction at the three-phase contact 



line. It is noted that this equation is identical to equation 
twenty- two of [33] for U = 0 and <tcos(0°) = <r sa - (r s i 
(Young's equation). Left hand side (LHS) of Equation 7 
is the out-of-equilibrium interfacial energy which is the 
driving force of capillary flow. Right hand side (RHS) of 
Equation 7 only includes dissipation of the free energy 
due to the contact line friction. De Ruijter et al. [30] 
showed that the corresponding dissipation function 
(TZi) is: 



qi 2 

2 



(8) 



In the next section, the wedge film viscous dissipation 
is calculated and added to Equation 8 to form the total 
dissipation function from which the total drag force is 
calculated. The total drag force is then equated to the 
LHS of Equation 7 to form the complete equation of the 
three-phase contact line motion. 

Hydrodynamic theory 

To calculate the wedge film viscous dissipation (TZ W ), 
Navier- Stokes equation of motion is solved in the wedge 
film region. From Figure 4 for the film thickness (H) 
much smaller than the radial distance p (H « p) and for 
capillary number Ca « 1, lubrication theory is used: 



dp d ( du\ 
dx = dzV n dz) 



(9) 



where p is the pressure and u is the velocity distribution 
at distance x inside the wedge film. For no stress bound- 
ary condition at the free fluid-air interface and no slip 
boundary condition at the solid surface, solution to 
Equation 9 gives: 



contact line region 



u(x^_ 


_H(p) 


Ae V 


P 


IV 




Figure 4 Schematic of a spreading droplet. 
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dp\ n n 



n \dX 



(10) 



c(cos(9 0 -cos(9) = (U 



+rj b F{cP)K 



2n + l 



l + n\l 



1-n J 6 n 



where rj n is replaced by its expression in Equation 1. 
The average cross-sectional fluid velocity in the wedge 

H 

film (u = | udz/H) is equal to the three-phase contact 



line velocity (u = U). This results in: 

2^ + 1 / / 7.\h+^ 



n + 1 



U 1- 1 



HJ 



(11) 



The viscous dissipation in the wedge film can be 
obtained as follows [5]: 



7T 



w 



K az 



r-x m n 

i j 

o L o 



dp 



(12) 



where r is the shear stress (= rj n d uld z), and x m is the 
cutoff length similar to slip length in HDT [27,28]. 
Without consideration of x m) dissipation of energy at the 
wedge film grows infinitely close to the three-phase con- 
tact line. For a thin wedge film, solution to Equation 12 
gives the wedge film viscous dissipation function: 



1-n 



(13) 



Dynamic contact angle 

Combining Equations 8 and 13 gives the total dissipation 
function [30]: 



U" 



1-n 
(14) 



Taking derivative of the total dissipation function with 
respect to contact line velocity (3 [TZi + TZ w ]/d U) re- 
sults in the total drag force [5]: 

/drag = 01 + „ b F(<p)K(^y(^J 1 [S-»- X J-»] W 

(15) 



Finally, equating Equation 15 with LHS of Equation 7 
gives: 



(16) 

It is noted that for n = 1 (Newtonian fluid), the integral 
of Equation 12 results in logarithm \n(r/x m ). I n this case 
the final form of Equation 16 is similar to De Ruijter's 
model [30] (<r(cos 6° - cos 6) = + 6n0(6)U\n{r/a)) 
where 0 = sin 3 (9/2 - 3 cos 6 + cos 3 6 and a is the cutoff 
length in De Ruijter's model). 

In Equation 16, the base radius (r) is in millimeter 
length scale while the cutoff length (x m ) is in nanometer 
length scale. Thus, r » x m , and consequently r 2 ~ n » xj^ n 
for n ranging from 0.04 to 0.92 (see Table 1). Also, for a 
sessile droplet of spherical geometry (see Figure 2), the 
base radius is geometrically related to the dynamic con- 
tact angle: 



3V 



sin 3 6 



n 2-3cos#+ cos 3 0 



(17) 



where V is the volume of the droplet. Neglecting x]^ n 
and substituting r with Equation 17 gives: 



cr(cos(9 0 - cos#) 
i:U + rj b F((P)K 



1-n 



2n + l 
n 

3V 



sin 3 6 



1 + n\ V tt 2-3 cos# + cos 3 <9 



1-n 

~3~ 



U n 



(18) 



Equation 18 shows the dynamic contact angle (6) as a 
function of contact line velocity (U), solid-liquid mo- 
lecular interactions (0, and non-Newtonian viscosity {n } 
K). Finally, substituting U with drldt = {drldO) x (dO/dt) 
the following equation can be obtained for the time evo- 
lution of the dynamic contact angle: 



a( cosa- cosa 0 ) = ( 



3V (1 + cosa) b 
n (2 + 3 cosa- cos 3 a) 4 

/ 



da 
dt 



2n + l\ n (\ 



1-n 



3V 



jt 2 + 3 cosa- cos 3 a 



V 



3V (1 + cosaf 
(2 + 3 cosa- cos 3 a) 4 



fda/dt 
V Ji-a . 



(19) 



Radiom et al. Nanoscale Research Letters 2013, 8:282 
http://www.nanoscalereslett.eom/content/8/1/282 



Page 6 of 9 



in which the dynamic contact angle 6 = n - a. To com- 
pare with experimental data 6 is used. Equation 19 is an 
implicit ordinary differential equation, which cannot be 
solved analytically, and thus numerical solutions to this 
equation will be sought. 

Results and discussion 

The effective diameter of nanoparticles was equal to 260 
nm at the lowest solution concentration of 0.05 vol.%. 
At higher particle concentrations, the increased 
interparticle interactions result in larger clusters. This 
increases the possibility of clusters to deposit on the sur- 
face of solid and form a new hydrophilic surface. Due to 
their larger size, these clusters are less possible to de- 
posit on the three-phase contact line, and thus a hetero- 
geneous surface will form: within the wedge film and 
away from the three-phase contact line, deposition of 
Ti0 2 clusters results in a hydrophilic surface with higher 
surface energy (approximately 2.2 J/m 2 [34]) than the 
three-phase contact line where the bare borosilicate 
glass is present (approximately 0.11 J/m 2 [35]). The 
higher surface energy inside the droplet shrinks the 
wetted area by increasing the equilibrium contact angle 
(denser solutions are more hydrophilic inside than 
outside). As a result, solid-liquid interfacial tension 
increases which on the other hand enhances the equilib- 
rium contact angle [21]. Surface tension of these solu- 
tions decreases with particle concentration that is in 
accordance with Gibb's adsorption isotherm. The shear 
thinning viscosity of the solutions is due to strong 
interparticle interaction of the nanoparticle clusters 
[19,23,36]. Other nanofluids such as ethylene glycol- 
based ZnO nanofluid [23] and CuO nanofluid [37] also 
exhibited shear thinning viscosity at low shear rates. 

Equation 19 suggests that the contact line friction dis- 
sipation (first term on the RHS of Equation 19) and the 
wedge film viscous dissipation (second term on the RHS 



of Equation 19) can occur at different time scales [38]. 
The time dependence of these dissipations has been 
shown by our experimental data: Figure 5 shows the ex- 
perimental three-phase contact line velocity (U = drldt) 
plotted versus o cos 0, where the base radius r is calcu- 
lated from the experimental dynamic contact angle 6 
using Equation 17. Figure 5 shows a linear trend that is 
in accordance with the contact line friction dissipation 
and a nonlinear trend (see inset of Figure 5) that is in 
accordance with the wedge film viscous dissipation. This 
suggests that at the start of capillary flow, the contact 
line friction is the dominant dissipative mechanism. As 
capillary flow slows down, the wedge film viscous dissi- 
pation becomes more dominant. This corresponds to 
the solution's higher viscosity at lower shear rates (see 
Figure 3). Transition to wedge film viscous dominant re- 
gime occurs earlier in dilute solutions; for example, 
Figure 6 shows that for 0.05% concentration the viscous 
forces start to dominate at time scales around 4 to 8 s 
while for 2% concentration at time scales around 25 
to 32 s. 

Figure 6 shows the dynamic contact angle of Ti0 2 -DI 
water nanofluids at various nanoparticle volume concen- 
trations ranging from 0.05% to 2%. Due to limitation in 
camera frame per second speed (30 fps), the onset of 
pendant droplet touching the surface of solid cannot be 
determined accurately. Hence, the time axis in Figure 6 
was shifted to where all of the captured images were 
readable to the FTA200 software. From Figure 6, it is ob- 
vious that for higher nanoparticle concentrations, the 
contact angles are higher. Figure 6 also shows that the 
spreading of these nanofluids starts from a primary re- 
gion where the contact angle changes rapidly followed 
by a region where the contact angle changes more grad- 
ually (note that in a very short period of time (less than 
300 ms), the contact angle evolves from 180° at point of 
contact to angles that are readable to our software and 
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Figure 6 Dynamic contact angle of Ti0 2 -Dl water solutions. 



are plotted in Figure 6 at the shifted zero time). In the 
primary region, the contact line friction dissipation 
predominates the wedge film viscous dissipation causing 
fast reduction in the contact angle; then the wedge film 
viscous dissipation controls the droplet spreading [31]. 

Using Equation 19, ( is obtained for the best fit of the- 
ory to experimental data that gives the least squared 
error. Figure 7 shows a reasonable comparison between 
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Figure 7 Dynamic contact angle of Ti0 2 -Dl water nanofluid, 
comparison of experiment and theory. 



experimental data and theory. The error is calculated 
from the following equation and is reported in Table 2: 



error 



\ 



E[(^-c periment )/ e ] : 



N 



(20) 



where N is the number of experimental points, Qf is the 
fitted dynamic contact angle, ^ x P eriment i s the experimen- 
tal dynamic contact angle, and e = 1 ° is the experimental 
error associated with the standard deviation in the meas- 
urement of the contact angles. A close examination of 
table three in [31] and table four in [38] reveals that the 
agreement between experiment and theory in our case is 
reasonable considering the complexity of the solution. 

Table 2 shows values of ( for various nanoparticle vol- 
ume concentrations. From solution concentration of 
0.05% to 0.5% C only changes by 5%; however, it drops 
rapidly for denser solutions. It is possible that the rela- 
tive higher hydrophobicity at the three-phase contact 
line for denser solutions lowers the affinity of surface 
molecules to water molecules, thereby lowering the fric- 
tion. At dense concentrations, the presence of large 
amount of nanoparticles in the wedge film varies the 
flow field structure. Without nanoparticles, it has been 
stated that there are two flow patterns in the wedge film: 
rolling and lubricating patterns [5]. Nanoparticles in the 
wedge film can change these flow patterns and result in 
more complex flow structures. As a result of these 
interparticle interactions, dissipation is more pro- 
nounced in the wedge film. Equation 19 gives better re- 
sults at lower nanoparticle concentrations since complex 
interparticle interactions are less frequent in dilute solu- 
tions (see Table 2). Other sources of disagreement be- 
tween experiment and theory can be local variations in 
the concentration of the nanoparticles in the nanofluid 
[21], pinning of the contact line, and variations in solid- 
liquid interfacial tension (a si ) [18,21]. It is not possible to 
model all these effects in theory, and only simple models 
which can accommodate some of these effects can be de- 
veloped. Also shown in Table 2 are the theoretical equi- 
librium contact angles, theory > wmcn are in reasonable 

Table 2 Coefficient of contact line friction C, theoretical 
equilibrium contact angle 6? heory , and error of comparison 
between theory and experiment 

Nanoparticle concentration 

2% 
1% 
0.5% 
0.1% 
0.05% 



C [Pa-s] 


e° 

^theory 


Error 


32 


52.1 


1.1 


99 


48.2 


1 


464 


46.4 


0.65 


483 


45.3 


0.54 


486 


44.8 


0.34 
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agreement with the experimental equilibrium contact an- 
gles, 6>° exp (see Table 1). 

Conclusions 

Due to a wide range of industrial applications, studying 
capillary flow of liquids laden with metallic and metal 
oxide nanoparticles is important. Metal oxide Ti0 2 
nanoparticles are especially interesting in enhanced heat 
removal applications. Agglomeration of nanoparticles re- 
sults in clusters that have larger effective diameter than 
the actual particle size. These clusters can deposit on the 
surface of solid substrates and form a heterogeneous 
surface condition inside the droplet away from the 
three-phase contact line that increases the equilibrium 
contact angle. Dynamic contact angle of metal oxide 
Ti0 2 nanoparticles dispersed in DI water revealed two 
stages of spreading: rapid reduction in contact angle co- 
incides with contact line friction dissipation governed by 
MKT while gradual reduction in contact angle coincide 
with wedge film viscous dissipation governed by HDT. 
Non-Newtonian viscosity of the solution is incorporated 
in HDT model to give reasonable comparison with ex- 
perimental data. Nanoparticles in the wedge film change 
lubricating and rolling flow patterns and result in com- 
plex flow field structures. Including all physical aspects 
of such complex flow in theory is not feasible at the 
current stage. Simple theoretical equations can only give 
reasonable comparisons with experiment. 
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